Semifiexible Polymers in a Random Environment 
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We present using simple scaling arguments and one step replica symmetry breaking a theory 
for the localization of semifiexible polymers in a quenched random environment. In contrast to 
completely flexible polymers, localization of semifiexible polymers depends not only on the details 
of the disorder but also on the ease with which polymers can bend. The interplay of these two 
effects can lead to the derealization of a localized polymer with an increase in either the disorder 
density or the stiffness. Our theory provides a general criterion for the derealization of polymers 
with varying degrees of flexibility and allows us to propose a phase diagram for the highly folded 
(localized) states of semifiexible polymers as a function of the disorder strength and chain rigidity. 
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PACS numbers: 36.20.-r Macromolecules and polymer molecules - 05. 40. -a Fluctuation phenomena, random 
processes, noise and Brownian motion - 75.10.Nr Spin-glass and other random models - 71.55.Jv Disordered 
structures; amorphous and glassy solids 
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I. INTRODUCTION 

In problems concerning flexible polymers in disordered 
media, disorder seems to play a more important role than 
the conformational properties of the polymers [1-12]. The 
effect is usually manifested in the size of the polymer, 
which in a simple but rather unrealistic case of Gaussian 
chains in a random potential, is found to be independent 
of the chain length L. The radius of gyration in three 
dimensions simply scales as R ~ b/(A/b 3 ), where & is a 
typical cut off length and plays the (undetermined) role 
of the segment size, and A is the disorder strength at 
a given system volume. The problem, however, is aca- 
demic mainly on two accounts: firstly, the effects due to 
excluded volume interactions are completely ignored, and 
secondly the effects of finite extensibility are not captured 
by the infinitely flexible Gaussian chains. It is conceiv- 
able that the inclusion of these effects, though nontrivial, 
can modify our present understanding of the problem. 

An effort in this direction has recently been made by 
accounting for the effects of excluded volume interactions 
[l2^ . The results show that a chain undergoes a der- 
ealization transition as the strength of the self-avoiding 
interaction is increased. Since localization is usually con- 
sidered as filling up of the entropic traps — created by 
randomness — with polymeric materials, it is under- 
standable that a self-avoiding chain can fill these holes 
up to a density close to pb 3 ~ 0(1) (p is the monomer 
density) before delocalizing and moving to the next hole. 
For Gaussian chains, on the contrary, derealization does 
not occur simply because they have an infinite stored 
length and no limit up to which the holes can be filled. 

In this paper we address the second of the above men- 
tioned problems, i.e., semifiexible polymers in a random 
environment. Since semifiexible polymers are finitely ex- 
tensible, the theory rectifies the infinitely flexible Gaus- 
sian chain model and serves as a useful description for 
localization of chains with varying degrees of flexibility. 
In contrast to Gaussian chains, localization of semifiex- 
ible polymers depends rather strongly on the conforma- 



tional properties of the polymer, i.e., the ease with which 
a chain can bend and eventually adopt highly packed 
folded states. While a polymer localizes at a certain fixed 
strength of the disorder, it can get delocalized with an 
increase in the stiffness which can be brought about by 
decrease in the temperature. A chain of a given fixed 
stiffness, on the other hand, can get delocalized with the 
increase in the disorder strength. The interplay of the 
disorder strength and bending energy is instrumental in 
bringing about localization of semifiexible polymers. 

The following section reviews the salient features of the 
wormlike chain (WLC) model; it provides a general cri- 
terion for localization of a semifiexible chain in random 
environment using simple scaling arguments. Section III 
summarizes the analytical difficulties involved in treat- 
ing the WLC model and discusses the Weiner model for 
stiff chains. Section IV presents a rigorous analysis of 
semifiexible chains in random media using the one step 
replica symmetry breaking calculations. The results are 
summarized in Section V. The details of the one step 
replica symmetry calculations are included in Appendix 
A. 



II. MODEL AND SCALING 

The configurational statistics of a semifiexible polymer 
is described by the wormlike chain model [15lll6| . In this 
model, a chain is considered as a continuous curve of total 
length L parameterized by an arc length variable s. The 
variation of the unit tangent vector u(s) = dr(s)/ds at 
different points along the curve with respect to the arc 
length variable determines the energy cost of bending [15- 
24]. To the lowest scalar order in du(s)/ds, the bending 
energy is given by 



i p ( L fdu 

2 i dS U 



(1) 



where [3 is {k^T) 1 and l p is the persistence length of the 
polymer. It is a measure of the distance over which bond 
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vectors are correlated, and hence determines the effective 
stiffness of the chain. In general, the correlation between 
the unit tangent vectors at two different points along the 
chain backbone decays exponentially, and is given by 



(u( S )-u( S '))=exp(-| S - S '|/«- 



(2) 



The persistence length is related to the bending energy 
e by l p = [3e. If such a chain interacts with TV randomly 
distributed obstacles, the Hamiltonian is of the form 



where r(s) is the position vector of a segment at distance 
s from one end of the chain; and is the position vector 
of the ith obstacle. W is the potential that describes 
the interaction between the polymer and the obstacles, 
the form of which is yet to be specified. If the obstacle 
density p(r) is defined by 
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P(r) = ^<5(r-r,), 



(4) 



the potential can simply be written as 

JV 



i=l 



drW[r(s) — r]p(r). 



(5) 



At the simplest level of description, the obstacle density 
can be assumed to have a Gaussian distribution [2J, |3j : 



/V(r.] = -I>(4/*^W^ 



where po is the mean obstacle density given by 



Po 



1 

V 



drp(r). 



(6) 



(7) 



To describe the effect of the disorder, one needs to per- 
form an average over the obstacle density. Since the av- 
erage is quenched, the estimation of the free energy relies 
on a difficult task of computing the average of the loga- 
rithm of the partition function Z . The problem is usu- 
ally circumvented by introducing n independent replicas 
of the initial system such that 



\nZ = lim(Z™ - l)/n. 



(8) 



The replicated partition function can easily be averaged 
over the Gaussian random potential given by Eq. (6) to 
yield 

(Z n ) = %[r( S )]exp[-/??U (9) 

Jr(0)=0 



where the effective Hamiltonian is given by 



m n = f E / ds 

a = l J° 



( <9u a (s) 



po 
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V ds 



E / ds ds'U[T a (s)-r (s')}(lO) 

a ,/3=l J ° J ° 



where 



U[r a (s)-r (} (s')]= dvW[v a {s)-v]W[T {s')-v] (11) 



If the interaction between the polymer and the obstacles 
is assumed to be short ranged, i.e., W[r] cx <5(r), then the 
potential can be expressed in terms of the natural scales, 
i.e., U(r) = glpS(r). In writing the last expression for 
U (r) , we have introduced a dimcnsionless coupling con- 
stant <?, which in general contains microscopic details. 
The natural (short) lengthscale defined in the problem is 
the persistence length l p itself. Therefore, we can express 
the scaling of the potential in terms of the persistence 
length; this makes the second term of Eq. (10) dimen- 
sionless, as should be the case 2]. Having substituted 
this form of the potential, the effective Hamiltonian is 
given by 



IE 

a = l 



ds 



du a (s) 



- E / ds 



ds 

L 



ds'6[r a (s)~r p (s%(12) 



where A = glpPo and has the dimensions of length. For- 
mally the second term in the above expression can be 
interpreted as a random potential which acts along the 
contour of the chain; it is a starting point of many studies 
on polymers in a disordered environment (see e.g. |13|). 
In the present study, the use of an array of obstacles in- 
stead of a random potential introduces a natural length 
scale in terms of the mean distance between the obsta- 

1 /3 

cles, lobs — 1/Po • It turns out to be useful since this 
length scale can directly be compared with the persis- 
tence length in various physical situations. 

Let us proceed by providing a rough estimate of the dif- 
ferent energies involved using simple scaling arguments. 
This can be done by dimensional analysis of the above 
Hamiltonian. In the case of annealed average there is no 
coupling between various replicas, that is a = 0, and the 
replica trick is not necessary. The free energy in three 
dimensions scales as 



am-h^L- a— 



(13) 



The first term accounts for the bending energy penalty 
(per k#T) for a chain of length L. It is typically defined 
as the ratio of the bending force constant (3e (= l p ) to 
the radius of curvature R. Note that the first term in 
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Eq. (|12|) corresponds to the bending energy per thermal 
energy of a semiflexible chain. The derivative of the tan- 
gent vector u(s) is nothing but the local curvature. The 
second term corresponds to excluded volume type inter- 
actions between the chain segments, which by virtue of 
the disorder average, are effectively attractive. In the an- 
nealed case the problem is similar to a semiflexible chain 
with attractive interactions. 

In the case of quenched disorder, on the other hand, 
there is a coupling between various replicas and the free 
energy scales as 



(14) 



where the second term in Eq. (|14fl describes a typical 
energy barrier of the quenched random potential; it is 
usually estimated as the root of the variance using sim- 
ple dimensional analysis 0, El ■ It is to be noted that 
A = A | In V |, where V is the volume of system. The 
issue of finite system volume is important for the case 
of quenched disorder since in an infinite random medium 
uenched and annealed disorders lead to the same results 
ITU . In what follows it is assumed that the system 
volume is much larger than the polymer dimensions i? 3 . 
The minimum of quenched free energy estimate yields 
the following result for the size of the chain, 



,4 



(15) 



which, as in the limit of totally flexible chains, does not 
depend on the chain length itself but is completely de- 
termined by the disorder and the persistence length. Al- 
though the size turns out to be similar to that of a Gaus- 
sian chain, its dependence on l p instead of b suggests a 
natural condition on the dimensionless ratio l p /A. That 
this indeed is the case can be seen from the following ar- 
guments: when the chain is reasonably stiff l p ss L and 
i?i oc can at most be of the order of L, which immediately 
implies that l p /L < A/l p ; in the limit of sufficiently flex- 



< R 



gauss 



ible chains, it is expected that i?i oc 
which implies {lp/L) 1 / 2 < A/l p . 

In addition to the upper cut-off on the size Ri oc , there 
is a lower cut-off, which is important for a flexible chain 
and suggests that i?i oc can not be less than the persis- 
tence length l p . R\ oc > l p therefore implies A/l p < 1. 
For a flexible chain, i.e., L/l p > 1, the latter constraint 
together with the constraint l p /L < A/l p can be summa- 
rized as 



1/2 



A 



(16) 



The second inequality implies that l p < (gpo) -1 ^ 3 , which 
suggests that as the chain become stiffcr it localizes at 
progressively lower density of the disorder po- In other 
words, since Z b s ~ Pq the chain is localized as long as 
its persistence length l p is less than the mean separation 



between the obstacles l b s - It is to be noted that the 
size of the chain has a simple interpretation in terms of 
the obstacle density. This can be seen by substituting 
A oc IpPo into Eq. l|15f) to produce i?i oc oc l p (l bs/lp) 3 - 
Then the criterion for localization, as discussed above, is 
simply given by the comparison between the persistence 
length and the mean distance between the obstacles. 

As discussed in the Introduction, a Gaussian chain un- 
dergoes a continuous decrease in its size with an increase 
in the obstacle density since the chain can be packed 
without limit into the cavities formed by the obstacles. 
In the presence of excluded volume interactions, however, 
a Gaussian chain can fill this cavity only up to a point 
before delocalizing and moving on to the next cavity; at 
high obstacle density the chain eventually localizes by 
forming pearl necklaces (as discussed in Ref. 12]). The 
present theory shows that a semiflexible chain without 
excluded volume interactions can dclocalize when the ob- 
stacle density and persistence length are related by Eq. 
(|16fl . This is in marked constrast to the infinitely flexi- 
ble Gaussian chain which can delocalize only when aided 
by the excluded volume interaction. In what follows we 
shall study the localization-delocalization transition of a 
semiflexible chain in the absence of excluded volume in- 
teractions, paying special attention to the flexible limit 
(L/lp > 1). We will show that even in the flexible limit, 
new physics arises because the persistence length defines 
a lower cut-off length scale which is absent in the Gaus- 
sian chain. We will also compare our results in this limit 
to earlier results for the Gaussian chain. 

Eq. I|16[l provides a general criterion for the derealiza- 
tion of a polymer with varying degree of flexibility. The 
next section addresses some of the difficulties involved in 
treating the wormlike chain model and presents a simpli- 
fied model to treat semiflexibility. 



III. GAUSSIAN STIFF CHAINS 

The wormlike chain model, though a proper descrip- 
tion of semiflexible polymers, suffers from its analyti- 
cal intractability in a large number of applications due 
to the constraint of finite extensibility, i.e., |u(s) 2 | = 
1, Vs G [0, L]. In many of the practical problems, there- 
fore, it is customary to relax the local constraint of unit 
tangent vector with the global constraint (|u(s) 2 |) = 1 
[17-24], producing a rather simple model in which the 
chain Hamiltonian has a connectivity term along with 
the bending energy term: 



dsu(s) 



ds 



du(s) 
ds 



(17) 



The mean square end-to-end distance obtained using the 
above Hamiltonian is identical with that of the worm- 
like chain model when 77 = 3/2l p and e = 3l p /2 [17-24]. 
The correlation between the unit tangent vectors, when 
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calculated from the above Hamiltonian, is of the form 



<u( S )-u( S ')>=exp(-| S - S '|/M: 



(18) 



where = 2Z p /3. The comparison of the above Equation 
with Eq. (2) suggests that the persistence length for this 
approximate model for stiff chains is smaller by a factor 
of 2/3 [23, l2ll l2~ij . This is not a serious concern in the 
context of the present problem since the Gaussian model 
contains the essence of semifiexiblity, i.e., finite extensi- 
bility, while remaining analytically tractable. Moreover, 
the difference between l v and lo becomes less important 
in higher dimensions |2l| . 

As before, we begin by writing down the free energy for 
this simplified model using dimensional analysis or Flory- 
Imry-Ma arguments which includes the Wiener stretch- 
ing term: 



IpL 



IpL 

W 




1/2 



(19) 



As expected, when the first term is ignored, and the free 
energy is minimized with respect to R, one obtains i?i oc = 
lp/A. However, the Wiener-term allows for more entropy, 
and it is natural to ask for perturbation around i?i oc . We 
look for solutions of the form 



R = Ri oc (l + 5), 



(20) 



where S is a small correction. When Eq. (16) is mini- 
mized with respect to R, and Eq. I|20|) is substituted in 
the resulting expression, 5 turns out to be of the form 



6 oc 



i+(i/y 2 (A/y 



(21) 



In the flexible limit of l p < L, Eq. Ijl6(l suggests A/l p < 1 
and Eq. (|2*U|l reduces to 



R^R 



loc 




(22) 



Although these scaling arguments appear to be reason- 
able and in accord with the physical intuition, the ap- 
propriate conditions and justifications remain to be re- 
covered from a more rigorous analysis. The next sec- 
tion presents a systematic approach to this problem us- 
ing Feynman variational method in replica space. Since 
for delocalized phases the replica symmetry breaking is 
no longer essential, the idea is to see if criteria for lo- 
calization and derealization can be derived from stan- 
dard replica techniques using one step replica symmetry 
breaking. 



IV. VARIATIONAL METHOD IN REPLICA 
SPACE 

A semiflexible polymer in presence of the randomly 
distributed obstacles confined in a harmonic potential is 



of the form 

PH. = '± \ ds 



n i L , fdr(s)\ 2 e f L , fd 2 r(s) 
1.11 'as 



ds 



ds 2 



N r L 



+ J2f daW[T(a)-Ti] + £ f dsr(s) 2 , (23) 
i= i Jo 1 Jo 

where u is the strength of the harmonic potential and 
defines the system size; the importance of this term will 
be discussed shortly. 

Following the same steps as used to obtain Eq. I|12|l . 
the quench average over the obstacle density produces 
the following effective Hamiltonian: 



IE / * 



dr a (s) 
ds 



d 2 r a ( s)\ | fj, 
ds 2 



fW dsv a {s) 2 

n r-L r-L 

V / ds \ ds'S[v a (s)-r b (s')}. (24) 



It is to be noted that carrying out an average over the dis- 
order density produces an effective attractive interaction 
of strength A between different replicas. To estimate the 
average free energy we employ the Feynman variational 
method in replica space [2, Uli , which is of the form 

n (F) = (H n h n ) hn - ~ J U n a=1 V[v a (s)]e-^, (25) 

where h n is the Gaussian trial Hamiltonian with varia- 
tional function er(s — s') characterizing the strength of 
the harmonic potential: 



V - f L 
0K = ~22 / ds 

0=1 J° 



IE . * 



a=l 



dr a {s) 
ds 

d 2 v a {s) 

ds 2 



a=l Ji) 

- ds ds'a ab (s - s')r a (s) ■ r b (s). (26) 

a,b=l J ° JO 

Defining the matrix of correlations by 

(r a ( S ) • r 6 ( S ')) = I ]T e- iqis ~ S,) Ga b (q), (27) 

and substituting into Eqs. l|5I|). (|2*5|l and the free 

energy simplifies to 

^j— = Const +7f-J^2^2 ( eq4 + 7 1<1 2 + V-)Gaa{q) 



a—1 q 



5 



2n ^ 

a, 6=1 



2n ^— ' /n 



(i.i: 



(2tt) 3 

dk 



(2tt) 3 



x e 



• S« ! ( G -(9)+ G "(9)- 2G ^ e "" !X ) 



(28) 



where delta function of Eq. I|24fl has been exponentiated 
and averaged over the Gaussian random potential to pro- 
duce the last two terms in the above equation. G ab (q) is 
the propagator with respect to the trial Hamiltonian h n : 

Gab(q) = [{eq 4 + + - M?)]" 1 - (29) 

&ab(q), in general, is an n x n matrix, which in the case 
of replica symmetric solution can be represented by a 
simple variational ansatz: 

<J a b{q) = -S ab [a' d S q . + a d (l-S q . Q )}+a(l-S ab )5 q . , (30) 

where the diagonal variational parameters oV's repre- 
senting the q — modes have been separated from 
all other modes (represented by <7<j's) to avoid working 
with infinite number of variational parameters. The off- 
diagonal terms have been considered to depend only on 
the q = mode, as is conventionally done Trans- 
lational invariance of the variational function requires 
J2 b a ab(q = 0) = 0) resulting in a condition a c { + a = 
that leaves only two of the variational parameters inde- 
pendent. 

It has been argued by Cates and Ball that in the ther- 
modynamic limit of infinite system size the annealed and 
quenched average should be the same, and an infinite 
Gaussian chain should always be collapsed Q| . That this 
indeed is the case has been shown by Goldshmidt, where 
the u — ► limit corresponds to an infinite system size, 
and the replica symmetric solution for a quenched ran- 
dom potential recovers the results of the annealed disor- 
der. On the contrary, when /i is taken to be finite, the 
replica symmetric solution is found to be an inadequate 
description of the problem, especially at high disorder 
strength where replica symmetry needs to be broken • 
The importance of weak replica symmetry breaking was 
recognized by one of us earlier |10| , where a special field 
theoretic description of polymers in strong disorder was 
used. In general, for problems concerning random po- 
tentials with short range correlations, one step replica 
symmetry breaking has been found to be a sufficient de- 
scription [25|. In what follows, therefore, we employ one 
step replica symmetry breaking to estimate the size of a 
scmiflexible polymer in a quenched random distribution 
of obstacles. 

From Eq. (27), the mean square end-to-end distance 
of a polymer in terms of the propagator G(q) is given by 

(R 2 ) = lim - V / ^(G aa (q)+G bb (q)-2G ab (q)e-^ L ), 
two n ' / 27T 

a,b=l J 

(31) 



where in the limit of large L the summation over q has 
been replaced by an integral. Eqs. (29) and (30) when 
substituted into the above equation followed by integra- 
tion over q produces 



<R 2 ) = 



1 fa 1/2 c- a - 2 " 1/2L/£l/2 

1/2 1/2 ^ C 



(32) 



where a± = [1 ± (1 — Aea d /rj 2 ) 1 / 2 ]/2 and ad is the vari- 
ational parameter to be determined by minimizing the 
variational free energy within one step replica breaking 
scheme fillip. 

In replica symmetry breaking scheme, the nxn matrix 
G(q) in the limit of n — > is parameterized by a function 
g(q,u), where u is a continuous variable that is defined 
in the interval [0, 1]. The diagonal elements, on the other 
hand, are parameterized by g(q) such that 



lim n G «b(q) = g(q) - (g(q)) 



(33) 



a,b=l 



where (g(q)) = J Q dug(q,u). 

The idea is to compute the inverse of a hierarchical 



matrix H(q) = [(eq 4 



^)S a b ~ v a b{q)] such that 



G(q) = Hiq)- 1 . The matrix H(q), in the limit of n -> 0, 
is similarly parameterized by the diagonal element h(q) 
and off-diagonal element h(q, u) given by 

h(q,u) = eq 4 + rjq 2 + (7 d ' S q . Q + (l-S q . )(7 d -a(q,u), (34) 

In one step replica symmetry breaking, the symmetry is 
broken at a point u c , and the function is defined by 



cr(q,u) = (T 5 qfi ; u < u c , 
= <7i S q . ; u > u c . 



(35) 



Translational invariance of the variational function in the 
limit of n — * for the case of one step replica symmetry 
breaking requires a d ' + u c ao + (1 — u c )ai = 0. The vari- 
ational free energy, as given by Eq. (|28|) . in the limit of 
n — > is therefore of the form: 



0(F) 



L 



const 



3 f dq 4 
2 J ^ q 



w 2 + fi)g(q) 



lim of / P-TrHH(q)} 
n^o 2n J 2n 



— — I d.r 



4 / ,//; 



dk 



dx 



-k 2 / !f9(9)(l-e- s ' x ) 



dk 

J2mf 



(36) 



The details of the replica symmetry breaking (RSB) 
method are described in the paper by Mezard and Parisi 
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|25| . the appendix of which provides the explicit expres- 
sions for some of the terms introduced in the above equa- 
tion. One step replica symmetry breaking steps relevant 
to the present problem are furnished in Appendix A. The 
free energy as given by Eq. (A14), in the limit of large L 
and fi — > 0, is dominated by 



const- 



3a d go 



-A 1 / 2 



9a 



3/4 



(37) 



where g Q = (77 + 2a 1 J 2 e 1 / 2 )- 1 / 2 ^^ 1 / 2 and A = 
3A| ln^|/(27r) 3 / 2 . The minimization of the free energy 
with respect to a d produces the following algebraic equa- 
tion: 



A\ V + 2a 1 J 2 e^ 2 y-a d = 0. 



(38) 



That Eq. I|32|) is a valid description in both stiff and 
flexible limit can be seen as follows: In the absence of 
disorder (corresponding to a d — > 0), when Li] 1 ! 2 /e 1 / 2 3> 
1, (R 2 ) w 3L/tj(= 2Ll p ) (the flexible limit); when 

v X/2 L / e l/2 < 1; ^ R 2^ ^ L 2 ( the ro( j lim j t ) 



V. RESULTS AND DISCUSSION 



DELOCALIZED 




FIG. 1: Numerical solution of Eq. 1391 representing local- 
ization - derealization transition as the dimensionless vari- 
ational parameter cro and dimensionless disorder strength 
Aq = A/io are varied. 



A. Localization-Delocalization Transition 

When i] = 3/2i p and e = 3i p /2 are substituted into Eq. 
(38), it can be rewritten in the following dimensionless 
form: 



Ao 4 (l + 3 ( 7o 1 / 2 ) 7 -a = 0, 



(39) 



where Ao = A/io, °o = cr d lo 3 and Iq — 2l p /3. 

The range of applicability of one-step RSB solution 
can qualitatively be seen from the behavior of u c , which 
is defined in the interval [0, 1]. The form of u c , as given 

by Eq. (All), suggests that |ln//| < LA 1 / 2 /g^ 4 , where 

.90 = l 2 (l + 3<t q )- 1/2 (t~ 1/2 . Since |ln/x| ~ \lnV\, the 
above inequality implies that a chain localizes when the 
binding energy is greater than the translational energy. 
It is to be noted that this inequality arises as a natural 
condition on the validity of replica symmetry breaking 
solution. In the flexible limit of L/Iq > 1 and Iq — * 0, Eq. 
(39) suggests that Co ~ (A/io) 4 , and the above inequality 
yields (Z /£) 1/2 | hi^ 1 / 2 < A// , the form of which is 
same as suggested using simple scaling arguments after 
Eq. IJT5J. 

Eq. (39) can be solved numerically for ao as a function 
of Ao. The results are summarized in Fig. (1); the de- 
localized region corresponds to the case where one-step 
RSB does not provide physically meaningful solutions to 
the dimensionless variational parameter cto . As is evident 
from Fig. (1), the chain is localized as long as the condi- 
tion A/Zo < 1 is satisfied. The latter condition introduces 
an upper cut-off on the dimensionless ratio A/io; it sug- 
gests that in the flexible limit (0 < l p /L < 1), for a given 



persistence length, there exist a critical disorder strength 
A c (< io) above which the chain delocalizes. 

The Gaussian approximation for semiflexible chains 
forbids us to consider rigorously the localization of a rigid 
rod [17-21]. In what follows, we restrict ourselves to the 
limit of intermediate stiffness < l„/L < 1. The real 
stiff limit of L/l p < 1 requires a separate treatment. 



B. The size of a semiflexible polymer in disorder 

In the flexible limit of io — > 0, Eq. (39) can be solved to 
quadratic order in 0$. The resulting expression is given 
by 



1/2 



A4[1 + (1 + 12(1-A4)/7A4)V2] 



2(1 -AS) (40) 

When 77 = 3/2i p and e = 3i p /2 arc substituted into Eq. 
(31), the resulting expression can be written in the fol- 
lowing form: 



<R 2 ) 



3Z 2 



x (a 



(1 + 3^)1/2 al /2 



1/2 ^2o 1 _ /2 L/3; 



1 



1 



1/2 



a+ — a 
1/2 



1/2 



where a± = [1 ± (1 - 9<t ) 1/2 ]/2. In the flexible limit 
of L/lo > 1, to leading order in io, the above expression 
can be simplified to produce the following approximate 
expression for the mean square end-to-end distance: 



3i 



(R 2 ) = ^(l-exp(-a i/2 L/i )) 



1/ 



(42) 
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FIG. 2: Log-log plot of the normalized mean square end- 
to-end distance (R- 2 ) /3Ll p vs. the dimensionless disorder 
strength Ao = A/£q. The curve is the numerical solution of 
Eqs. £23 and JUJ for the ratio L/l = 100.0. 



In the limit of low disorder strength Ao — > 0, R « lo/A, 
which is same as Eq. (I15|) obtained using scaling argu- 
ments. When the disorder strength reaches its critical 
value corresponding to Ao « 1, R w ^o/A 2 - The lat- 
ter suggests that when the average separation between 
the obstacles is of the order of the persistence length 

— 1/3 

(^obs ~ Pq ~ lo), the disorder can explore the length 
scale over which bending energy is stored. Since any fur- 
ther increase in the disorder density amounts to high en- 
ergy penalty for bending, the polymer delocalizes above 
the critical density. The decrease in the size close to 
A/lo ~ 1 can therefore be attributed to the dominance 
of the bending mode at this length scale. The results are 
summarized quantitatively in Fig. (2), which is a numer- 
ical solution of Eqs. (|39J) and Ij41|l and shows a sharp 
decrease in the size of the polymer close to the critical 
disorder strength. These results are in contrast to to the 
conventional studies on the localization of a Gaussian 
chain, which is a fractal, and undergoes a monotonous 
decrease in the size with the increase in the disorder den- 
sity. In the limit of A/Iq < 1, to the first order correction 
in A/lo, the size is given by R w l 2 (l — A A /lf/)/A, which 
is similar to Eq. (|22Jl . The one step replica symmetry 
breaking, therefore, supports the scaling arguments pre- 
sented in Section III. 

It is to be noted that in the above analysis we have 
used Iq instead of l p as the persistence length. This is 
because ^o — 2i p /3 appears as the natural persistence 
length for Gaussian stiff chains, an issue that has been 
discussed in Sec. III. In the rest of our analysis l p will be 
used as the persistence length. 



C. The final state of the polymer 

When the polymer is able to find a suitable region 
among the obstacles, its conformation must contain a 
rich structure. Though the results of the previous sec- 
tion show a dramatic decrease in the chain size with the 
increase in the disorder strength, they say nothing about 
several different folded states a chain can adopt as a func- 
tion of the chain rigidity L/l p and dimensionless disorder 
strength A/l p . In what follows we propose a phase dia- 
gram for several possible compact states of the localized 
polymer by comparing their surface energies. 

Since an average over the disorder density produces an 
effective attractive interaction between chain segments, a 
chain can find itself in several compact states depending 
on its stiffness. For a fully flexible chain the compact 
state is expected to be a globule; a stiffer chain, on the 
other hand, has a tendenc y to wind around itself many 
times to form a toroid ptll27|. In presence of obstacles, 
the size of a circular toroid can be obtained by minimizing 
the free energy of the form 

r-^-iff •« <43 > 

where b is the monomer size and l p is the length over 
which monomers of size b are correlated. The first term 
accounts for the bending energy of a toroid 28]; R s is 
the radius of the torus tube, and R is the radius from the 
center of the hole to the center of the tube. Minimization 
of the free energy with respect to R yields R ~ l p /A 

and R s ~ A 1 / 2 L 1 / 2 b/l p , where the volume constraint, 
V = 2n 2 RR 2 = Lb 2 , has been used. As expected, the 
minor radius R s grows at the expense of the major radius 
R with the decrease in the persistence length implying 
that for a given disorder strength there exist a critical 
ratio of L C rit /lp above which the chain prefers to be in a 
globular state. Comparison of the surface energy penalty 
of a toroid St — "fRRs with respect to a sphere S g ~ jR 2 
suggests L crit /l p ~ (l p / A) 3 (l p /b) 2 , where 7 is the surface 
tension. As expected, for high disorder strength, A/l p , 
and, L/l p > L cr i t /l p , the chain prefers to be in a globular 
state. In the opposite limit the chain adopts a toroidal 
state. 

When the disorder strength is high, a chain can pack 
densely only if it chooses a mean direction that a given 
distribution of obstacles imposes. In other words, a high 
obstacle density favors anisotropy in structures, and it 
seems unlikely that under these conditions a polymer can 
adopt uniformly bent states like globules or toroids. It is 
physically intuitive, however, that if the polymer is able 
to localize for moderate persistence lengths, it can do so 
by producing many hairpins in a small available volume. 
The hairpins are the abrupt reversals in the direction of 
the chain from the mean ordering direction and result 
from the balance between the elastic penalty that favors 
gradual bends, and orientational field that favors rapid 
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DELOCALIZED 




GLOBULES 



persistence length; n max — L/X determines the maxi- 
mum number of hairpins that a chain of a fixed contour 
length L can sustain [2fj ISfT ] . By comparing the surface 
energy cost for hairpins, Sh — jl p R, with that of toroids 
one obtains L crlt /l p ~ (A/l p ) 3 (X/b) 2 (X/l p ) 2 . In writing 
the surface energy cost for a hairpin, we have accounted 
for the fact that the average length of the chain in di- 
rections parallel to the mean ordering direction scales as 
l p . It turns out then that for high disorder strength and 
L/lp > L CT n/lp a chain prefers to form many hairpins. 
In the opposite limit toroidal state is favored. The re- 
sults are summarized in the form of a schematic phase 
diagram in Fig. (3). It clearly shows the dominance of 
hairpins at moderate persistence length and high disorder 
strength; for small and large persistence lengths the chain 
conformation is dominated by uniformly bent shapes like 
globules and toroids respectively. 



FIG. 3: A schematic phase diagram based on the scaling anal- 
ysis depicting different folded states of a localized polymer as 
a function of the dimensionless disorder strength, A/l p , and 
chain rigidity, l p /L. The coexistence curves are the regions of 
comparable surface energies. At large l p /L for a given large 
A/lp the chain is delocalized in the static sense since under 
these conditions (Eq. I|16|l i the mean separation between the 
obstacles is less than the persistence length and the chain 
cannot find an appropriate place in the disordered medium to 
sufficiently bend; it can do so only at the cost of very high 
bending energy. 

reversal in the alignment of the chain in the preferred 
direction. The energy costs for forming a hairpin is given 
by 

PUh = {u h l P ) 1/2 , (44) 

where Uh is the coupling constant of a (mean-field) Maier- 
Saupe potential between the aligned chain segments |29| . 
Since Uh depends on the microscpic details, a simple di- 
mensional analysis of the above equation shows that Uh 
can at most be of the order of 1/6. The free energy for 
this case is given by 

^=^W/ 2 - f^J L, (45) 

where we have assumed that the number of hairpins rih 
must be resolved in a self-consistent way, i.e., = L/R. 
By minimizing the the free energy with respect to R, we 
find for the final state of the localized chain 

where A represents the length of a bend given by A = 
y/lp/uh, and is of the order of less than or equal to the 



VI. CONCLUSIONS 



We interpret our results to suggest that the interplay 
between the two length scales — the persistence length 
of the chain and the average separation between the ob- 
stacles — serves to localize a semiflexible polymer in a 
random distribution of obstacles. As long as the dis- 
order density is low, the mean separation between the 
obstacles is larger than the persistence length, and the 
chain is localized; but as soon as the disorder explores 
the length scales of the order of the persistence length 
the chain delocalizes since the energy penalty for bend- 
ing exceeds the energy gain in binding. On the contrary, 
when the density of obstacles is kept fixed, derealiza- 
tion can be brought about by varying the temperature. 
While the ease of localization of a semiflexible chain can 
be enhanced by increase in the temperature (correspond- 
ing to the decrease in the persistence length), a localized 
flexible chain can be delocalized by lowering the temper- 
ature. Once the polymer is localized it adopts a highly 
compact folded state by forming hairpins for moderate 
persistence length and high disorder strength; globules 
and toroids are favored for small and large persistence 
lengths respectively. 

Although the idea was to present a general theory for 
localization of semiflexible polymers, there is one respect 
in which the present formalism could be seen to have lim- 
itations, and that is the use of the Gaussian approxima- 
tion for stiff chains, which fails to adequately describe the 
rigid rod limit. Nevertheless, given the range of systems 
that can be described by flexible chain models, this the- 
ory holds out the possibility of addressing several realis- 
tic problems, for instance, the immobilization of biopoly- 
mers on surfaces with random interactions. 
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APPENDIX A: ONE-STEP REPLICA 
SYMMETRY BREAKING 

The definitions of some of the terms introduced in Eq. 
is as follows: 



lim i-Trln[H(q)} = ln(%) - (%))) 
n— >0 n 



h(q,0) 



h(q) - <%)) 



Uu h(q)-(h(q))~[h](q,u) 
— ^ In 



h(q) 



~9(q) 



1 



h(q) - (h(q)) 



h(q,0) 



h(q) - (%)) 
[h](q,u) 



o ^h(q)-(h(q))-[h](q,u) 



(Al) 



(A2) 



i(q,u) = -~ 



l 



h(q) - (h(q)) 
" dv 



[h](q,u) 



u(h(q)-(h(q))-[h](q,u)) 
[h](q, u) 



o v 2 h(q)-(h(q))-[h](q,u) 
%,0) 



h(q) - (h(q)) 



(A3) 



For the case of one step replica symmetry breaking the 
calculation is straightforward, and the explicit expres- 
sions for the above equations are given by 



h(q) - (h(q)) 

m 

g(q,u < u c ) 
g(q,u > u c ) 
with 



eq 4 + r 1 q 2 + n + a d (l-S qfi ), (A4) 
| __ | Si 

/' a* 2 V 

(1 - S q>0 ) 



= S q ,0 



(eq 4 + nq 2 + a d ) 



= Sq,0 



Si 



CO 



u c fi(n + 'E 1 ) /i 2 



(A5) 



(A6) 



[h](q, u) = ; u < u c 

= -Si S qfi ; u> u c 



(A7) 



where Si = u c (ai — uq). Having substituted the above 
equations into Eq. (|3l)|l followed by integrations with 



respect to the variables q and k, the expression for the 
free energy simplifies to 



0(7) 



const 
3 

+ 2l 



3(m - o- d )go 



2g (ea d )^ 



x lln 
AL 



u, 
Si 



\i /i + Si 
1 



2(2^)3/2 5q 3 / 2 V7o 

- g(u > u c y 3/2 ) - u c (g(u < u c )~ 3/2 
- > u c )- 3 / 2 )) , (A8) 

where go — (rj + 2o- 1 J 2 e L / 2 )~ 1 t 2 o- d ~ 1 / 2 ] other terms intro- 
duced in the above expression have the following defini- 
tions: 



2/1 1 

g(u >u c ) = 1 + — — — 

L 9o \M + s i a d 

2/1 1 

g(u < u c ) = 1 + 



Lg \fj, + Si a d 
Si 



U c Ai(/i + Si) 



(A9) 



(A10) 



1 



fi( x ) = 1_ ^72 172 ( 



1/2 -a^^/'Lx/e 1 ^ 



a+ e 



aj_ — a_ 

1/2 -a^^Lx/e^'s 



(All) 



where a ± = [1±(1 — ^e(T d /rj 2 ) 1 / 2 ]/2. The free energy can 
be minimized with respect to Si and u c . In the limit of 
large L and fj, — ► 0, the solution of the resulting equations 
yield 



Si = AVV 7/4 ; 



(A12) 
(A13) 



where A = 3A| ln/x|/(27r) 3 / 2 . The substitution of above 
equations into Eq. (|36[) implies the expression for the 
free energy to 



(7) 3(p - a d )g 



1/2 -3/4 



X -3/2 r r l 

A.g L ' 



dx( gi (x)-^ 2 - 1) 
6|ln/x| Jo 

_f^2 /X-l/a 7/4 _ -is 

2|ln M | 1 50 d j ' 



(A14) 



where only terms of 0(1) with respect to expansion in 
1/L have been retained. 



10 



[1] 


A. Baumgartner and M. Muthukumar, in Polymeric Sys- 


[17] 


R. A. Harris and J. E. Hearst, J. Chem. Phys. 44, 2595 




tems, vol. XCIV of Advances in Chemical Physics, ed. by 




(1966). 




I. Progogine and S.A. Rice, p.625 (Wiley, N.Y., 1996) 


[18] 


K. F. Freed, J. Chem. Phys. 54, 1453 (1971); M. G. 


[2] 


S.F. Edwards and M. Muthukumar, J. Chem. Phys. 89, 




Bawendi and K. F. Freed, ibid. 83, 2491 (1985). 


2435 (1988). 


[19] 


K. F. Freed, Adv. Chem. Phys. 22, 1 (1972). 


[3] 


A. Baumgartner and M. Muthukumar, J. Chem. Phys. 


[20] 


B.-Y. Ha and D. Thirumalai, J. Chem. Phys. 103, 9408 


87, 3082 (1987). 




(1995). 


[4] 


M.E. Cates and R.C. Ball, J. Phys. (Paris) 49, 2009 


[21] 


M. Otto, J. Eckert and T. A. Vilgis, Macromol Theor 


(1988). 




Simul 3, 543 (1994). 


[5] 


J. Machta, Phys. Rev. A 40, 1720 (1989). 


[22] 


S. M. Bhattacharji and M. Muthukumar, J. Chem. Phys. 


[6] 


J. Machta and T. R. Kirkpatrick, Phys. Rev. A 41, 5345 




86, 411 (1987). 




(1990). 


[23] 


J. B. Lagowski, J. Noolandi and B. Nickel, J. Chem. 


[7] 


T. Nattermann and W. Renz, Phys. Rev. A 40, 4675 




Phys. 95, 1266 (1991). 




(1989). 


[24] 


R. G. Winkler, P. Reineker and L. Harnau, J. Chem Phys. 


[8] 


S.P. Obukhov, Phys. Rev. A 42, 2015 (1990). 




101, 8119 (1994). 


[9] 


S. Stepanow, J. Phys. A 25, 6187 (1992). 


[25] 


M. Mezard and G. Parisi, J. Phys. I, 1, 809 (1991). 


[10] 


P. Haronska and TA. Vilgis, J. Chem. Phys. 101, 3104 


[26] 


A. Y. Grosberg, Biophysics, 31, 1471 (1991). 




(1994). 


[27] 


J. Ubbink and T. Odijk, Biophys. J., 68, 54 (1995). 


[11] 


Y.Y.Goldschmidt, Phys. Rev. E 61, 1729 (2000). 


[28] 


J. N. Bright and D. R. M. Williams, Europhys. Lett., 45, 


[12] 


Y.Y.Goldschmidt and Y. Shiferaw, Eur. Phys. J. B 32, 




321 (1999). 




87 (2003). 


[29] 


M. Warner, J. M. F. Gunn and A. B. Baumgartner, J. 


[13] 


G. Migliorini, V.G. Rostiashvili, T.A. Vilgis, Eur. Phys. J. 




Phys. A, 18, 3007 (1985). 




B 33, 66 (2003). 


[30] 


D. R. M. Williams and M. Warner, J. Phys. France 51, 


[14] 


S. Stepanow, J. Phys. (Paris) I 2, 273 (1992). 




317 (1990). 


[15] 


H. Yamakawa, Modern Theory of Polymer Solutions, 


[31] 


D. R. M. Williams and A. Halperin, Macromolecules 26, 




Harper and Row, New York, 1971. 




2025 (1993). 


[16] N. Saito, K. Takahashi and Y. Yunoki, J. Phys. Soc. 








Japan, bf 22, 219 (1967). 







